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ABSTRACT 


The low-speed bluff body studies conducted on aeroelastic 
models were extended by redesign of the installation in the West 
Coast Research Corporation's subsonic wind tunnel at the U. S. 
Naval Postgraduate School, Monterey, California. Improvements 
included stiffer support base, redesign of the spring constraint, 
and the use of eddy-current type displacement gages. [In addition, 
the model's lateral response data were monitored, conditioned and 
recorded on one-inch magnetic tape in a voltage analog form. 

The analog data were converted to digital form at a sampling 
rate of 500 numbers per second for a period of eight seconds per 
case in order to permit subsequent digital analysis in the Control 
Data Corporation 1604 computer. The results of this digital 
computation provided a numerical estimate of the spectral 
distribution, which gives an insight into the character of the 
random response of the aeroelastic structure due to the 


aerodynamically induced unsteady air load. 
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ШІ; Introduction. 

The use of aeroelastic type models has been widespread and is 
well known for estimating flutter boundaries on complete aircraft. 
Early airplane flutter models were primarily concerned with matching 
vibration frequencies, mode shapes, and mass scaling without 
concern for compressibility or viscosity effects. The advent of 
supersonic aircraft brought on an impetus to scale Mach number, 
hence the design and testing techniques on aeroelastic models 
changed in accord with the usage of transonic and supersonic 
wind tunnels [1]. The fairly recent importance of ground wind 
loads on missiles has brought emphasis upon bluff-body type 
aeroelastic models with similar scaling factors, except that 
Reynolds number is substituted for Mach number. Hence, a revision 
occurs in testing missile elastic models since pressurized low- 
speed wind tunnels are important in order to provide proper 
Reynolds number matching. 

In addition to studying specific three-dimensional bluff bodies, 
a need exists for determining and evaluating both the forcing 
functions and structural response characteristics of idealized 
single degree of freedom bluff body aeroelastic models with cross- 
sections as a test parameter. Consistent with this type of work, 

a program started in 1965 by Breckon and Hite [10] has been extended 
and the test installation and technique has been upgraded in order 
to provide a better degree of insight into bluff body model response. 
Engineering-type improvements included a stiffer support base, 


redesign of the spring constraint, and the use of eddy-current type 
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displacement gages. These gages operate as variable inductance 
instruments where inductance changes are a function of air gap. 

Specifically, with regard to the antenna project the horizontal 
members, termed boom boards, have exhibited undesirable wind induced 
oscillations. Ideally the boom board structure could be visualized 
as a constant cross-section beam having some unknown end fixity 
(approximately a simple support) with an evenly spaced elastic load 
over its length and acted upon by an unknown random aerodynamic 
forcing function resulting from steady wind or gusts. 

The boom boards have cross sections of circular cylinders, 
squares, and rectangles. They all have the common feature of being 
bluff bodies, which means that the flow does not progress smoothly 
over the body, but separates from the surface in an unsteady manner 
which gives rise to the generation of fluctuating air loads. RMS 
amplitude of the response increases with increased wind velocity. 
The lightly damped structure acts as a band-pass filter to the 
random forcing function [2]. Usually the structure responds only 
to the forcing function centered about the resonant frequency of 
the structure. The principle response is essentially at the body's 
natural frequency. 

The response of the model differs from that of the usual 
flutter model since bluff body response is associated with random 
type signals, while airplane flutter models involve harmonic type 
outputs. The analysis of random signals can be approached by pure 
analog signal processing methods where the random signal is 


recorded on magnetic tape and then played back through a harmonic 
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analyzer. Drawbacks to this technique have usually been that high 
quality filter circuits are needed when operating in the two cycle 
per second band width range; the lower limit on frequencies is 
between 20 and 30 cycles per second, and only power spectral density 
results are obtained. The auto-correlation function has not yet 
been obtained very satisfactorily by pure analog means. 

Another technique used in handling random data is the hybrid 
approach; this method was utilized in this survey. The data were 
recorded in analog voltage form and then converted to digital form 
utilizing a sampling rate consistent with the frequency content of 
the signals. A sample rate of 500 numbers per second allows 
spectral distributions up to 250 cycles per second, while a 
digital record of eight seconds in duration provides a lower bound 
on spectra of about one cycle per second. The conversion to 
digital form utilized an analog/digital converter in conjunction 
with a Control Data Corporation (C.D.C.) 160 computer. Once in 
digital form, the binary tapes were processed through a C,D.C. 1604 
computer. The output results were estimates of auto-correlation 
and power spectral densities. 

The low-speed wind tunnel used for these investigations was 
readjusted and calibrated during the summer of 1965 and the flow 
traits were improved from those reported by Breckon and Hite [10]. 
Main changes included removing the flow swirl and inclination in 
the test section by means of corner turning vane adjustments. A 
five orifice yaw head was designed and constructed for this purpose 


and allowed measurements of flow angularity to +0.02 degrees. 
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Turbulence measurements were made in February, 1966 using a 
0.0002 inch diameter tungsten hot wire probe. In addition to the 
usual RMS measurements, the turbulence measurements were processed 
through the analog es digital data handling system in order to 
provide a measure of the scale of turbulence and the spectral 
distribution of random flow eddies. Although the addition of fine 
mesh screens in the tunnel settling section are desirable, the 
tunnel in its present condition is quite usable for measurements 
of a character like those reported herein. 

The author worked very closely with Lieutenant W. P. Mortenson 
throughout the entire antenna vibration survey; however, separate 
papers are being presented. The purpose of this paper is to 
specifically describe the aeroelastic model used to simulate the 
antenna boom boards,the data handling system utilized in data 
processing for the antenna vibration project, and to present the 
computer program for determining the power spectral dénsity, thereby 
defining the response of the aero elastic wind tunnel model to the 
random aerodynamic forcing function. The detailed analysis and 
explanation of the results provided by this system and recommenda- 
tions for "in the field fixes" for the antenna project can be 
founden | IIIF 

The investigation was conducted by the author in conjunction 
with Lieutenant W. P. Mortenson during the period October 1965 
through May 1966 at the U. S. Naval Postgraduate School, Monterey, 
California. My sincere appreciation is extended to the faculty 


and staff of the Aeronautical Engineering Department for their 
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assistance. The author is deeply indebted to Professor L. V. Schmidt 
for his invaluable encouragement and assistance throughout the project. 
23 Description of the Model and Equipment. 

In order to simulate the boom boards in the wind tunnel, it was 
necessary to design an aeroelastic model. The characteristics of 
a dynamic model are quite stringent. Ideally the model should be 
light weight, mounted on frictionless bearings, able to reproduce 
desired frequencies, and sufficiently stiff so as not to introduce 
undesirable structural oscillation. The general approach used was 
to approximate one half of the prototype by hinging the model at 
one end with the other end free, with controllable stiffness. Since 
boom boards are quite long and their end fixity conditions are 
close to simple supports, the first lateral vibration mode was 
assumed dominant. The aeroelastic model can thus be restricted to 
one degree of freedom in the lateral direction. 

Since the intent was to observe many different configurations, 
interchange of the model shape was a necessity. For simplicity, it 
was decided to make the center of the model the strength member. 

An extruded aluminum "I" beam was selected primarily for its light 
weight and ease of machining. This aluminum beam extended the | 
entire length of the model. 

The outside configuration was formed by shaping balsa wood to 
the desired dimensions. The balsa covering was constructed in two 
symmetric pieces and contoured to fit snuggly on the "I" beam 
(Figure 1). End plates were screwed into the balsa wood at the 


top and bottom; these end plates were mounted to the "I'' beam by 
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a bracket. The end plates were adjustable to ensure that the model 
was centered between top and bottom of the tunnel with only a small 
end gap remaining. 

The wind tunnel floor and ceiling were fitted with aluminum 
plates which had the center section recessed slightly to accept 
the end plates. This allowed the model to span the entire height 
of the tunnel and thus reduced the flow around the ends. 

When scaling a dynamic model for wind tunnel testing, the 
variables which must be considered are: fluid density and viscosity, 
fluid velocity, vibration frequency, and inertia effects (stifíness, 
mode shapes, mass distribution, etc.). In low velocity flow the 
fluid may be considered incompressible, i.e., Mach number effects 
may be neglected. 

The first lateral vibration mode has been assumed dominant, 
thus permitting the aeroelastic model to be restricted to one 
degree of freedom in the lateral direction. This allows the 
simplification that model inertia effects need not be scaled with 
respect to prototype. The remaining variables combine into the 
dimensionless parameters: Reynolds number (R=vd/y ) and 
Strouhal number (S= fd/v). 

Matching of Reynolds numbers between prototype and model 
assures that the aerodynamic forcing function, which causes the 
vibration, correspond. Matching of Reynolds number will effectively 
occur as far as dynamic model response is concerned when the air 


flow about the model is supercritical L6 719 Matching of Strouhal 
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number or dimensionless frequency assures dynamic response 
correspondence between prototype and model [8]. 

To illustrate these relations a typical example follows. 
Assume that the prototype bluff body has a characteristic depth 
of 12 inches, a fundamental resonance frequency of five cps, and 
responds to sea level velocities of 25, 50 and 75 miles per hour. 


Case Velocity Strouhal Reynolds 


mph {ps Number Number 
A 2 SSO? 0.136 2 .34x10? 
B 5009379 0.068 4.66x10 
C 75 110.0 0.046 7.00x10 


To determine the appropriate aeroelastic model to represent 
the prototype, enter Figure 2 with Reynolds number (R) and model 
characteristic depth (d) to determine the dynamic pressure (q). 
With this information enter Fígure 3 with q, d, Strouhal number (S) 
to determine the required resonance frequency of the model. Typical 


models which would match the above prototype condition are as follows: 


Case d q f 
inches psf cps 
A 50 2.5 30 
4.0 14.4 45 
B 6.0 29,2 20 
5.0 36.0 29 
4.0 ວໆາ 45 
C 6.0 57.5 20 


Instead of making an exact match on both Reynolds number and 
Strouhal number, certain classes of bluff body cross sections may 
be expected to not be as severely dependent upon Reynolds number 
and hence only require Strouhal number matching. A clue in this 


direction is provided by Figures 5 and 8 of the work by Delany 
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and Sorensen [5] where they present steady state drag coefficient 
as a function of Reynolds number for several cross-section shapes. 
The circular cross section is well known to have a transition with 
respect to drag coefficient at a Reynolds number value of about 
9 Ξ T 
3 x 10 while square and rectangular shapes do not exhibit this 
5 5 6 

trait at Reynolds numbers well in excess of 1 x 10. 

If we presume that the unsteady air loads behave in a like 
manner to the steady state drag coefficients, then the necessity 
of matching Reynolds number exactly in an aeroelastic model with 
a square cross section loses importance. On this assumption, the 
sample values considered would indicate the necessity to select 
model operating conditions such that Strouhal numbers of 0.05 to 
0.15 corresponded to the model resonance condition for Reynolds 


numbers of approximately 4.0 x 107. This condition could be 
matched by the model in any of the following configurations: 

d=6.0 inches, q= 18.5 psf, £ =12 to 38 CPS, OF 

d=5.0 inches, q= 21.7 psf, £ =16 to πο ος от 

d=4.0 inches, q= 42.0 psf, Ё =29 to 80 cps 

The top of the model, "Т" beam, was attached to a rigid support 

by means of a flexure plate (Figure 4). This flexure plate minimized 
frictional damping, while allowing motion primarily in the lateral 
direction. Further, it provided some drag restraint. An additional 
drag restraint was provided near the bottom of the model by means 
of a stranded stainless steel cable with turnbuckle arrangement 


for tensioning. 
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The bottom of the model was supported in a manner which provided 
the copmility of varying the model's natural frequency. A heat 
treated one-half inch diameter steel drill rod was attached to the 
aluminum "I'' beam by two bolts machined to accept ball bearings. 

The bolts were adjustable to allow centering of the drill rod, and 

the ball bearings provided a relatively frictionless pin joint 

(Figure 5). The drill rod was cantilevered to a movable steel 

support block, precision drilled to accept the spring rod, and 
provided with six set screws which fixed the bar in position. The 
movable steel support block was fitted to a 16 inch precision 

machined wedge shaped track. Eight lag bolts secured the sliding 
block in any position (Figures 6 A,B,C). The track was bolted to 

a one-inch thick steel plate supported by four welded four-inch 
diameter steel pipes which formed legs in a manner similar to a 
dynafocal suspension principle (Figure 7). This rebuilt model 

support base increased the rigidity of the model's reference base 

over that used in [10]. The improvement in the lateral stiffness 
constraint provided a means for more rapid adjustment of the spring 
bar length while assuring that the base of the constraint was attached 
to the support base. Figure 8 illustrates the complete specifications 
of the model. 

The wind tunnel utilized was the 42 inch by 60 inch subsonic 
tunnel located at the U. S. Naval Postgraduate School, Monterey, 
California. This tunnel is a closed section, vertical return, 

200 mile per hour installation built by West-Coast Research 


Corporation of Los Angeles, California. Power is supplied by two 
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150 horsepower electric motors driving two independent constant 

speed variable pitch fans. However, for ease of operation only one 
fan was utilized; this fan provided adequate power for the dynamic 
pressure required. In addition,a survey of the dynamic pressure 

on tunnel centerline was completed. In order to check the turbulence 
level of the tunnel, a hot-wire survey was performed. The results of 
these surveys are presented in Appendix A. The aeroelastic model 

was mounted vertically in the test section as shown in Figure 9 

along with an illustration of the entire data handling system. 

The primary instrument for response measurement was the Bently 
displacement measuring device. This device provides an instanta- 
neous measurement of the lateral displacement of the bottom of 
the beam. 

The Bently sensor unit, a small pancake shaped coil, is 
excited at a one megacycle carrier frequency by a transistorized 
amplifier unit and power supply combination. The rapidly 
fluctuating magnetic field at the face of the coil allows sensing 
of the air gap change between the coil face and a flat non-magnetic 
conducting surface due to eddy current effect. The frequency limit 
of the Bently is 10,000 cycles per second. This is well above the 
125 cycles per second anticipated during the antenna survey. In 
practice this device proved quite reliable with a consistent distance 
resolution of about 20 micro-inches. Figure 6B shows the Bently 
sensor unit mounted to the model support base. Therefore, the gage 
measured relative displacement between the support floor and the 


model. Appendix B provides a calibration for the Bently inductance 
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displacement measuring device. For test purposes, the air gap was 
allowed to vary in the linear range from 0.02 to 0.08 inches, where 
the output sensitivity was approximately 70 volts/inch. 

The secondary source of Feces frequency response data was the 
Gulton type A321 accelerometer. This piezoelectric type device was 
used with an Endevco amplifier. The accelerometer has a flat 
response curve to above 3,000 cycles per second. Results of Gulton 
accelerometer calibrations are presented in Appendix B. It was 
decided to use the Bently as the primary source of data due to its 
reliability and ability to measure relative displacement between 
fixed reference and the base of the model instead of the inertial 
reference frame of the accelerometer (Figure 6C). 

Figure 10 presents a block diagram of the entire data 
acquisition system. Due to the complexity of the wiring and the 
desire to have an orderly system, BNC fittings were used throughout 
in conjunction with shielded coaxial cable. The coaxial cable 
decreased the possibility of picking up any outside disturbances 
such as 60 cycle from power supplies, ground loops, and 
flourescent lights. 

All data were monitored and controlled through a central switch 
panel in a manner similar to a telephone patch panel in order to 
maintain a systematic control of data flow. Voltage signals from 
the Bently and accelerometer were monitored and recorded using a 
Ballantine model 320A voltmeter in order to obtain the true RMS 
reading of the randomly varying signal. This meter uses a diode 


squaring circuit so as to give an accurate measure of the RMS 
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voltage, and can handle situations where peak values occur up to 
ten times greater than the RMS. 

The Vidar, integrating digital voltmeter, reads both frequency 
and voltage. It was utilized in the system to insure that the 
initial setting on the Bently, with the model stationary, was in 
the center of the linear range of the response characteristics of 
the Bently. Its secondary function was to calibrate the center 
frequency of the tape recorders frequency modulating (FM) 
data channels. 

The tape recorder/reproduce unit used was the Ampex FR-100. 
This provided 14 channels for the recording of data. Both frequency 
modulated and direct record amplifiers were available to record 
data. With a recorder/reproduce speed of seven and one-half inches 
per second, the frequency response for the frequency modulated 
amplifier was from zero to 1250 cycles per second. The input to the 
recorder was limited to one volt RMS or Е 1.40 volt peak to peak. 
The voltage signals from the Bently and accelerometer were recorded 
through the frequency modulated amplifier. Voice announcements 
were recorded through a direct response amplifier. These 
announcements properly labeled and identified the data. The system 
recorded the data on one inch mylar instrument grade tape so that 
the voltage signals which represent an analog of the model response 
could be analyzed at a later date. 

Тһе Fektronic 545A cathode ray oscilloscope is a single beam 
device which behaves in a dual beam manner using the dual trace 


amplifier channel, and it was used to monitor various signals. 
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The dual trace provided the means for simultaneously observing the 
signals from both the Bently and the accelerometer. Since the input 
to the Ampex tape recorder was limited to one volt RMS, observing 
the input determined whether E а ол was necessary to keep the 
signal within limits. Finally, the scope could be connected to She 
Ampex reproduce amplifier (output) to insure that the data was being 
properly recorded and reproduced. 

3% Test Procedures. 

The entire data acquisition system with a few exceptions 
consisted of tube elements. Therefore, prior to obtaining data, the 
entire system was allowed a proper warmup time in order to insure 
that all thermal transients had settled out. The Ampex tape 
recorder/reproduce amplifiers were calibrated in accordance with 
the operation manual. 

Calibrations were recorded on the magnetic tape through the 
frequency modulated amplifiers on the exact channels that the 
Bently and acceleration voltage signals were to be recorded. This 
consisted of zero volts and plus and minus 1.4 volts (DC). This 
calibration could be utilized later during the data analysis phase. 

To determine the system damping and to calibrate the spring rod 
length as a function of the natural frequency of the model, static 
tests were performed with the wind tunnel inoperative. The model 
base was mechanically displaced from the centered equilibrium 
position and sharply released. The model's decaying response of 


both displacement and acceleration was recorded onto magnetic tape. 
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The vibrational response of the model was obtained by a series 
of test runs. Each was run at a fixed natural frequency set by 
positioning the calibrated spring rod to the proper length while 
varying the tunnel dynamic pressure. Here again, the data from 
the Bently and accelerometer were recorded on one-inch magnetic 
tape and the RMS values from the Ballantine voltmeter were read and 
recorded in the data logs. Visual observations, using the Tektronic 
oscilloscope, were made of the response wave form to insure that 
the amplitude did*not exceed record amplifier limitations, and 
further to check that data was being accurately recorded 
and reproduced. 

Voice announcements on a separate channel preceded and followed 
the calibration check, model frequency response check and the 
dynamic data runs. This insured that the data was properly labeled 
and would definitely be retrievable. 

The data from the one-inch magnetic tape was next played back 
through an analog representation of an averaging circuit and plotted 
to observe the continuity of the curves, and further to look for any 
violent or sporatic portions. As expected from the oscilloscope 
observations, the curves were smooth and well-behaved, although 
complex and random in nature, Further, it appeared that the data 
was a stationery random process (this was subsequently confirmed 
by power spectral analysis.) Since the data was well-behaved, any 
portion of the data could be analyzed without fear of large errors. 


The model frequency response check and portions of the data were 
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plotted on strip chart Brush oscillograph to obtain frequency. 
Samples of this data may be found in Appendix C. 
4. Data Processing System. 

At this point, approximately 450 data runs in analog form 
were recorded on one-inch magnetic tape. It was desirable to form 
the auto-correlation function and then the power spectral density 
in order to obtain the response traits of the various configurations 
to the random forcing function. It was estimated that it would take 
approximately three to five man-months to process each run by hand. 
Obviously it was necessary to utilize the Control Data Corporation 
(C.D.C.) 1604 digital computer. However, a major hurdle still 
existed, namely the conversion of the analog data into digital form, 
acceptable as computer input. 

The process of converting real analog data into computer 
digital form was accomplished in two distinct and separate steps. 

In the first step the analog data was digitized at discrete points 

and written in BIT form on one-half inch magnetic tape. The second 
step was to unpack the digital data and put it into a form useable 

in the C.D.C. 1604 as standard program input. 

Figure D-l is a block diagram of the analog to digital 
conversion system. The wind tunnel data from the one-inch magnetic 
tape was reproduced through the Electrical Engineering Department's 
Ampex CP-100, a solid state record/reproduce system. The digitizing 
computer was the CDC-160 located in the Electrical Engineering 
Department Computer Laboratory. The CDC-160 selects either the 


12 BIT analog to digital converter unit for input, or the CDC-163 
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tape unit for output. The computer has a 4096 cell 12 BIT memory 
length. The lower 4,000 cells were reserved for storage, and the 
top 96 cells were reserved for programming. The core storage was 
an obvious limiting factor. However, 4,000 data samples were 
assumed adequate. The basic concept of the operation was to take 
samples until the core storage was full (4,000 samples) and then 
write this information onto one-half inch magnetic tape. 

The procedure utilized a manually initiated pulse to trigger 
the analog to digital converter. The computer controls the inter- 
sample time delay and the analog to digital converter samples rapidly 
until 4,000 discrete samples have been taken and stored in the memory 
core. The computer then dumps the output onto one-half inch magnetic 
tape. The analog to digital converter then remains hung up until 
receipt of another manual pulse. Program Digitize, a CDC-160 tape 
program, is found in Table D-l. 

A sample rate of 500 discrete points per second allows spectral 
distributions up to 250 cycles per second, while a digital record of 
eight seconds duration provides a lower bound on spectra of 
approximately one cycle per second. This frequency range is more 
than adequate since the natural frequencies of aeroelastic models 
tested varied from 35 to 90 cycles per second. 

In order to identify each block of 4,000 data points on the 
tape, a 16 digit identifier was first written followed by the data 
points. These identifiers are extremely important since they are 
the only means by which the data may be retrieved at a later time. 


This 16 digit identifier is a manual input to Program Digitize. 
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The data in digital form may be recalled by subroutine 
Data. Subroutine Data, Table E-2, in turn uses subroutine Findit, 
Table E-3, which locates the identifier on the tape and subroutine 
Unpack, Table E-4, a machine language routine to unpack the samples 
so that one-to-one correlation exists between the original CDC-160 
listing and the unpacked CDC 1604 version. 

A short CDC-160 tape program Test, Table D-2, was used to 
adjust the bias voltage to minus five volts and the desired inter- 
sample delay to the proper time. Additional more detailed 
information of the digitizing system can be found in Lieutenant 
N. S. Barrett's thesis "Extracting Analog Signals From Noise Using 
The Digital Computer" [9]. 

In order to test the capabilities of the digitizing system, a 
test tape was prepared with a 250 cycle per second sine wave, square 
wave, triangular wave, 20 cycle per second band-limited white noise 
and finally white noise. The data was analyzed and presented 
graphically, Figures D2-6. The system worked successfully and 
provided very satisfactory results. 

The final step in the data handling system was the computer 
program, written in Fortran 63, which utilized the Tukey method 
for estimating the power spectral density [4]. 

The complete program may be found in Appendix E. The program `° 
possesses sufficient comments to make it self-explanatory. However, 


some additional comments on the input cards are provided in Figure E-l. 
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Basically, this program used the data recalling system 
previously described. The integer format) from the subroutines 
te to floating point format and converted into volts 
by a constant. In order to insure that no bias existed, the data 
were normalized to wipe out the DC term on an average basis. 

The auto-correlation function of the random response was 
computed, and subsequently used, to estimate the spectral distribution 
of the energy by applying a Fourier cosine transform to the auto- 
correlation function. For any given signal, aperiodic or random in 
nature, it is quite clear that the auto-correlation function and 
spectral density are quite meaningful, but only the auto-correlation 
function is directly measurable. The power spectral density may be 
defined as the energy distribution on an average basis in the 
frequency domain; that is, energy with respect to frequency. 
Appendix F gives a more detailed mathematical explanation of the 
auto-correlation function and power spectral density to show the 
inter-relationships. The final step in the program was to determine 
the energy contained in the power spectrum by application of the 
trapezoidal rule. Comparing this result with the total energy 
provides a rapid indication of the accuracy of the results obtained. 
A sample listing of the data output is given in Table E-5 and is 
graphically presented in Figure E-2. The data and graphs specifically 
define the frequency response of the aeroelastic model to the 
aerodynamic random forcing function. This information is imperative 
to the investigation of the aerodynamic phenomena and subsequent 


analysis described in [11]. 
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5% Summary. 

A controlled test environment for conducting wind-tunnel 
investigations of bluff body aeroelastic models has been established. 
The data acquisition and processing systems are precise, orderly, 
and operational. The conversion from an analog voltage record to a 
digital format has been accomplished with a high degree of accuracy. 
The computer program was utilized to analyze approximately 450 wind 
tunnel data runs, and presented the output in numeric and graphic 
form. These results were imperative to the success of the 
investigation and analysis of aeroelastic phenomena conducted on 


the various configurations as presented in [11]. 
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APPENDIX A 


WIND TUNNEL TEST SECTION CALIBRATION 


The test section pressure survey shown in Figure A-l was 
performed at the U. S. Naval Postgraduate School's low speed wind 
tunnel. The profiles shown were taken 22 inches aft of the test 
section leading edge, with reference dynamic pressure of 40 pounds 
per square foot. 

These wind tunnel flow traits were improved with respect to 
tunnel dynamic pressure and flow inclination. Although the 
importance of these changes upon bluff body separation may not be 
critical or discernible, it is good practice to have a well 
established and uniform test environment in the wind tunnel 


test section. 


Aconstant temperature hot wire was placed on tunnel centerline 
at the model test station in order to make clear tunnel turbulence 
measurements. The hot wire was 0.001 inch diameter tungsten, 
approximately 0.08 inches in length and was supported on the end 
of two needles projecting from a support tube. It was possible 
to make direct measurements of turbulence level relative to free 
stream velocity using a prototype constant temperature amplifier 
circuit recently developed at the U. S. Naval Postgraduate School. 
These turbulence measurements were read directly by the Ballantine 


True RMS voltmeter, and also recorded on magnetic tape for 


subsequent digital analysis using procedures described elsewhere. 
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The direct measurements of turbulence indicated a turbulence 
level of approximately 0.45%, but the single wire sensed 
fluctuations both in the vertical and axial direction. On the 
assumption of isotropic turbulence, it is possible to estimate 
the axial velocity turbulence level as approximately 0.32%. This 
value was obtained uniformly for a test section q range of 10 to 
40 psf. The above value of turbulence level is considered high 
for high calibre wind tunnel work in a low turbulence environment, 
but considering that no anti-turbulence screens were installed 
in the settling section, and that the test was concerned with 
bluff body separation, the value was quite acceptable. The 
installation of screens in the settling section should result in a 
much lower value of turbulence level. 

The results of the digital calculations are shown on Figure A-2 
with respect to both auto-correlation and power spectral density. 
The area of the mormalized auto-correlation curve is 0.102 seconds 
which may be interpreted as a measure of the scale of turbulence, 
that is, on the average, turbulent fluctuations are unrelated 
for time intervals beyond 0.102 seconds. The power spectral 
density curve shows no dominant energy peaks other than one at 
1/2 cps. It is quite possible that the 1/2 cps peak was caused 
by the vibration of the hot wire support. The absence of other 
peaks indicates an absence of noise due to fan blade rotation or 
an undamped flow separation from other sources in the tunnel 


circuit. 
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The basic conclusion reached from the hot wire measurements 
is that the test section circuit is functional in its present 
configuration and the addition of screens at a later date will 
make the facility very useful for fundamental low turbulence 
experiments. 

The tunnel airspeed and/or dynamic pressure was set using 
clear tunnel values. Іп making a final evaluation, 'solid blockage 
effects may be taken into account. An estimate of blockage for 


the case of a cylinder spanning the tunnel is given by the formula: 


where a = radius of the cylinder 


h = height 
E fractional velocity increase 
Ru M 5 a 
4 .0146 1.03 
Sa .0296 1.06 
8 .0584 1.12 
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APPENDIX B 


CALIBRATION OF BENTLY AND GULTON INSTRUMENTS 


The Bently sensor head was mounted in an aluminum test 
block. The procedure used was to obtain the gage voltage for 
a specific distance measured by a micrometer. The plot of 
distances versus output voltage is shown in Figure B-l. The 
slope of the linear region was 70 volts/inch. The Bently sensor 
unit was serial number 14184, and the distance detector unit 
model D-252 was serial number 5254. The ease of calibration and 
the consistent reliability made the Bently displacement measuring 
System an extremely useful instrument. 

The Gulton model A321, a piezoelectric accelerometer, is 
extremely sensitive and because of its high natural frequency 
it can be used through a wide frequency range. However, since 
it is an A.C. instrument it cannot be statically calibrated. 
Therefore, the accelerometer was calibrated against a previously 
calibrated Statham on the Calidyne shaker table. The Gulton 
accelerometer serial number 1558 was combined with an Endevco 
amplifier serial number 128. The calibration is shown on 


Figure B-2. 
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FIGURE B=2 


CALIBRATION OF GULTON ACCELEROMETER 
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APPENDIX C 


BRUSH OSCILLOGRAPH TRACES OF DATA 
AND FREQUENCY RESPONSES 


Typical brush recorder tracings of response data of the 
aeroelastic model to the random forcing function are :shown on 
Figure C-l. The lever arm and average frequency are shown on 
each trace. 

Typical brush recorder tracings of the static frequency 
response test are shown in Figure C-2. These are used to determine 


the system natural frequency (£5) and damping (B ). 
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FIGURE C-1 
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FIGURE C-2 


TYPICAL BRUSH RECORDER TRACINGS OF 
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APPENDIX D 


ANALOG TO DIGITAL CONVERSION SYSTEM 


The process of converting real analog data into computer 
digital form was accomplished by the system shown in Figure D-l. 
The system utilized the C.D.C. 160 computer with Program 
Digitize, Table D-1, to sample the analog data digitize and 
write the digital form on magnetic tape compatible with the 
C.D.C. 1604 computer. Program Test, Table D-2, is used to 
calibrate the system to insure accuracy. The data in digital 
form may be recalled by utilizing subroutine Data, Table E-2, 
which in turn utilizes subroutine Findit, Table E-3, which 
locates the pre-set identifier on the tape and subroutine 
Unpack, Table E-4, a machine language routine which unpacks 
the data. In order to test the capabilities of the digitizing 
System a test tape was prepared ,| The first 160 points of each 


is plotted in Figures D-2A through D-2E. 


Sil 


0101 
0603 
7064 
7500 
2410 
7500 
1401 
7600 
4176 
2076 
3465 
6134 
7500 
2400 
2074 
4160 
2067 
4161 
2075 
4162 
2073 
4163 
7500 
2111 
7303 
0000 
6102 
0155 
2246 
4076 
2073 
3473 
06155 
0401 
4075 
5455 
3454 
6153 
6061 
5476 
2066 
0701 
6501 
7056 
0760 
0000 


PTA 
ADNO3 
JP164 
EXFOO 
2410 
EXFOO 
1401 
INA 
STI76 
LDD 76 
SBD65 
NZF 34 
EXFOO 
2400 
LDD74 
STI60 
LDD67 
STT61 
LDD75 
STI62 
LDD73 
STI63 
EXF00 
2111 
OUTO3 
0000 
NZFO2 
0133 
LDF46 
STD76 
LDD75 
SBD73 
NZF55 
LDN01 
STD75 
AQD55 
SBD54 
NZF53 
ZJF61 
AQD76 
LDD66 
SBNOL 
NZBO1 
JPI56 
0760 
0000 
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TABLE D-1 


PROGRAM DIGITIZE 


New’ No? Na? Nr ft ft МУ 


Na Na? Na Nr rs 


Jump *to 
Initial 


Set enable 


A/D Ch. 1 


60 to 47 


Clear enable 


Form indent 


1011... 
Output From 
LWA + 1 


0137 A 
Сӣ бо 0115 
Set block No. = 1 


GØ TØ 0120 
GØ TO D127 


Jetas 188815 


Max blocks 760 
Total block count 


121 


122 
123 
124 
125 
126 
127 
130 
131 


132 


0005 
0003 
0185 
0134 
0135 
0136 
0100 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 


0000 
4077 
2200 
0137 
4076 
3072 
4031 
0701 
4065 
2070 
4074 
0401 
4075 
7077 
5475 
5455 
7057 
2074 
3471 


6103 
2074 
7701 
5474 
7057 
7500 
1111 
2055 


7702 


TARIE D= (Con't.) 


0005 
0003 
0133 
0134 
0135 


Gio EG = =. 


STD77 
LDCOO 
0137 

STD 76 
ADD 72 
STD31 
SBNOI 
STD65 
LDD 70 
STD74 
LDNO1 
STD/5 
JD 7 
AOD75 
AOD55 
JEN Sy 
LDD 74 
SBD71 


NZFO3 
LDD 74 
MLTO1 
A の D74 
JPI57 
EXF00 
ШІП 

LDD55 


HLTO2 


53 


lst word address 

2nd of header 

3rd of address 

LWA of header 

Addr of initial 

LWA of data 

Inter sample delay 

Spare ID (anything) 

Initial run no 

No runs wanted 

No samples/block 

No blocks/run 

Current run number 

Current block number 

Running Storage 
Address 

Jump control 

Enter initial 


Set FWA of data 
Compute LWA + 1 


Set LWA of data 
Load initial run no. 
Set I. R. N. 


Set l1l. B. BN. 

Exit initial 
Update block no. 
Update block count 
Go to 0003 


Check no. runs 
enough 


Current run no. 
Update run no. 


Go to 0003 


Display total 
block count 


TABLE D-1 (Con't.) 


133 0000 Run no. 10 
134 0000 Spare ID 
155 0000 Block ID 
136 0000 No. blocks/run ID 
1373 
) Data storage area 
7776 ) 
Note: M = Manual entry 


e 
lI 


Computed entry 
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-ມ 6. + UN + O 


7500 
2410 
7500 
1401 
7600 
4070 
7500 
2401 
73:08 
0072 
6102 
0070 
0400 
4071 
6414 
7700 
7500 
2410 
7500 
1401 
7600 
4143 
0300 
0300 
0300 
0300 
0300 
0300 
0300 
0300 
2042 
0701 
6501 
7044 


0043 
0022 
7701 


TABLE D-2 


PROGRAM TEST-160 


EXFOO 
2410 
EXFOO 
1401 
INA 
STD70 
EXFOO 
2401 
OUTO3 
0072 
NZF02 
0070 
LDN00 
STD71 
EJB14 
HLTOO 
EXFOO 
2410 
EXFOO 
1401 
INA 
STI43 
N@P 
N の P 
N の P 
N の P 
N の P 
N の P 
NOP 
N@P 
LDD42 
SBNOI 
ΝΖΒΟΙ 
JPI44 


Set enable 


Call A/D CH 1 


Output 
From 0071 to 0072 


Return to 000Z 
Set enable 


Call A/D CH 1 


Timing control word 


0043 


HLTO1 


S 


FIGURE Del 


DIAGRAM OF ANALOG TO DIGITAL CONVERSION SYSTEM 
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jJdBl 1831 1611910 КОНУ 
2180 LAdNI JO SLINIOd ALXIS GSYGNTH 1€81J 9838930 
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Jjdol 1531 61190 ИОН 
σι συ LIGNI dO SiNiOd ALAIS а-ЯСМПН 15414 94-454 
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APPENDIX E 


POWER SPECTRAL DENSITY PROGRAM 


The energy distribution on an average basis in the 
frequency domain is defined by the power spectral density. À 
computer program utilizing the Tukey method for estimating 
spectral distribution was written - Table E-1. Tables E-2 
through 4 are subroutines which locate and unpack the data 
desired. The program has sufficient comments to make it self- 
explanatory. Figure E-l is an illustration of the control card 
utilized in the program. A sample listing of the data output is 
given in Table E-5 and graphically presented in Figure E-2. The 
data and graphs specifically define the frequency response of the 
aeroelastic model to the random aerodynamic forcing function. 

It will be noted in Table E-5 that the first 100 values of 
input data are listed in order to give visual assurance that the 
numbers are valid. The listing of auto-correlation is normalized 
with respect to its value at t- 0 and is specified as AN(TAU) 
versus Tau(SEC). The value A(0) shown above the listing is the 
actual mean square value of the auto-correlation at time shift 
T= 0 seconds, and is in units of REDE. This value may be 
related to the tape recorder voltage signal by the scaling factor 
of 5.0 volts (digitized) = 1.4 volts (recorded). The power 
spectral density is listed as XN(FREQ) versus FREQ(CPS) and is 
normalized such that the area under the curve of XN versus FREQ 


is unity. 
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TABLE ЕЈ. 
POWER SPECTRAL DENSITY COMPUTER PROGRAM 


-COOP,,JESBERG BOX Ј,1/1/0/49/5/15/25/Е/45 =54/5=50/6=51,15,20000 
“ҒТМ,91»Е» 
PROGRAM POWSPEC 
C TUKEY METHOD FOR ESTIMATING POWER SPECTRAL DENSITIES 
DIMENSION P(4000)»A(501)0»X(501) »FREQ( 501)»TAU(501)5»Q0(3000) 
DIMENSION KDATA(4000 92) »sIDENT (200) 
DIMENSION ITITLE(12) 
DIMENSION DATE(160)»SKIP(160). 
COMMON KDATA 
EQUIVALENCE (KDATA»P) 
С FIRST STEP IS DATA CONTROL CARD--DATA READ IN BY TAPE OR CARD , 
C FINO IF DATA IS ON TAPE--ITAPE=0 FOR TAPE»=1 FOR CARD DATA 
WRITE (69103) 
1 READ (59101) ITAPEsIDELTAsDELTATsISTART 9 IMAX oM 
WRITE (69104) 
WRITE(6?101) ITAPE6IDELTA5DELTAT5ISTART 6 IMAX 6M 
1Q=ISTART-1 
IADJ=IMAX-IQ 
2 IFCITAPE - 1) 60595 
C START DATA, CARD READ IN--USE THIS BRANCH FOR TEST PURPOSES 
C IF CARD DATA REMOVE STATEMENT 11 AND REPLACE WITH A CONTINUE 
5 READ (5»102) (P(I)»I=1»IMAX) 
WRITE (65105) 
WRITE (6»106) (P(I)9»Iz19»100) 


GO TO 30 
€ START TAPE DATA READ IN--USE THIS BRANCH FOR PRODUCTION 
c NUMHDR = NOe OF HEADINGS 


6 READ 110»NUMHDR 
DO 7 N=1sNUMHDR 
7 READ 111»IDENT(N) 
WRITE (69114) 
WRITE (69111) (IDENT(N) sNs1 sNUMHDR ) 
C MAIN DO LOOP ON ENTIRE PROGRAM 
DO 500 N=1»NUMHDR 
IDENT =IDENT (N) 
KLIST=1 
MAX = IMA X 
CALL DATA(IDENT+MAX»+KLIST+KFLAG) 
IF(KFLAG) 85998 
8 WRITE (60112) 
WRITE (69113) Ne IDENT 
GO TO 500 
9 00 11 J=l»IMAX 
[51 
10 P(J)=KDATA(J»1) 
1 1 P(J) =P(J) /( -409。6 ) 
WRITE (69105) 
WRITE (69106) (P(I)sI=1+100) 
GO TO 30 
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30 
31 


32 


33 


34 


35 


40 
41 
42 


43 


44 


45 


46 
50 


TABLE E-1 (Con't.) 


NORMALIZE INPUT DATA TO WIPE OUT D-C TERM ON AVERAGE BASIS 
ASUM=0.0 

DO 32 I=lsIADJ»sIDELTA 

ASUM=ASUM+P (1) 

CONTINUE 

AZ 1S NUMBER OF DATA SAMPLES BEING COMPUTED 
AZsFLOATF(CIADJ-1)/CIDELTA))*1, 0 

ASUM=( ASUM) /( AZ) 

DO 34 Is1eiADJsiDELTA 

P(1)=P(1)-ASUM 

CONTINUE 

ASIS=0.0 

DO 35 I=19+160 

DATE(1)=P(1) 

SKIP(I)=ASIS 

А515=А515+1,0 

CONTINUE 

FIND MEAN SQUARE AZERO IsE。 AUTOCORRELATION AT TAU=ZERO 
ASUMz=000 

DO 42 I=1l»+IADJ+IDELTA 

ASUM=ASUM+P (1 )**2 a 

CONTINUE \ 

AZERO=ASUM/AZ 

FIND AUTO-CORRELATIONS A(L) FOR Lz1,M 
DO 46 L=1»M 

ASUM=0,0 

BSUM=0.0 

CSUM20, 0 

MZ=(L*IDELTA+]) 

DO 45 I=MZ, IADJ»IDELTA 

IZ=-I-tL*IDELTA) 

ASUMZASUM-*P(IZ) *P(I) 

BSUM=BSUM+P (IZ) 

CSUM=CSUM+P (1) 

CONTINUE 
AZ2(1.0)/FLOATF(CCIADJ-1)/(IDELTA))*1-L) 
AC(L)ZCAZ) *XLASUM)-(AZ**2 ) *BSUMKCSUM 
CONTINUE 

FIND XZERO---POWER SPECTRAL DENSITY AT FREQsZERO 
А50М=0,0 

MZ =M-1 

FM=FLOATF(M) 

CS1=COSF(3。14159/FM) 
SN1=SINF(3。14159/FM) 

CSL =CS1 

SNL=SNI 
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51 


G 
С 
52 
C 
53 
54 


55 
56 
57 
58 
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ТАВІЕ Е-1 (Con't.) 


DO 52 L=1»MZ 
AZ=HANNING FACTOR 

FIND COSINE SUM FROM TRIGe IDENTITIES 
AZ=1。0+CSL 
ASUM=ASUM+AZ#A(L) 
CSL1sCSL*CS1-SNL*SN1 
SNL1=SNL#FCSI+CSL#SN] 
CSL=CSL1 

SNL=SNL1 

CONTINUE 

DZ=0。5/FM 

XZERO=DZ#( ASUM+AZERO) 
FIND X(K)-—--POWER SPECTRUM AT Kzx=loM 
CSK=CS1 

SNK=SN1 

рО 59 К=1,М 

А50Мя0.0 

CSKL=CSK C 
SNKL=SNK 

CSL=CS1 

SNL=SN1 

DO 55 L=1»MZ 
А7х(1.04С5( )%С5КІ. 
ASUM=ASUM+AZ#A(L) 
CSL1=CSLACS1-SNLASNI 
SNL1-SNL*CS14CSL*SN1 
CSL=CSL 1 

SNL=SNL 1 
CSKL1zCSKL*CSK-SNKL* SNK 
SNKL1zSNKL*CSK4*CSKL SNK 
CSKL=CSKL] 

SNKL=SNKL 1 

CONTINUE 
IF(K-M)56657657 

02=1,0 

GO TO 58 

DZ=0。5 

GO TO 58 

DZ =DZ/FM 
X(K)=DZX(ASUM+AZERO ) 
CSK1=CSK*CS1=SNK*SN1 
SNK1=SNK*¥CS14CSK#SN1 
CSK=CSK1 

SNK=SNK 1 

CONT INUE 
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60 


61 


62 


63 


65 


66 


67 


68 


TABLE E-1 (Con't.) 


SHIFT INDICES--ZERO FREQ WILL START AT K=l AFTER SHIFT 
USE STORAGE SPACE OF P(I) DATA DURING SHIFT OPERATION 
NORMALIZE ANUTOCORRELATION 

NORMALIZE AUTOCORRELATIONWsReTs AZERO 

DO 61 1=1»M 

P(I+1) =X(1) 

P(I+1001)=A(1) 

CONTINUE 

X(1)=XZERO 

A(1)=160 

М2 =М+1 

DO 63 1=2»MZ 

X(1)=P(1) 

A(1)=P(I+1000)/(AZERO) 

CONTINUE 

APPLY TRAP. RULE TO FIND ENERGY CONTAINED IN POWER SPECTRUM 
FOR RANGE OF K=1»M+1---DEFINE ENERGY AS XENGY 

NOTE X(K) IS ENERGY (WeReTe UNIT CHANGE OF INDEX K 
ASUM=0.0 

DO 66 К=2,М 

ASUM=ASUM+X (K) 

CONTINUE 

XENGY=0。5*(X(1)+2。*ASUM+X(M+1 ) ) 

FIND FRACTION OF TOTAL ENERGY IN CALCULATED FREQ。 RANGE 
XFACT=XENGY/AZERO 

OBTAIN SPECTRAL DENSITY WeReT, FREQ(CPS)--NORMALIZE WeReT XENGY 
FREQ(1 ) =0。0 

ТАО(1)-0.0 

MZ=M+1 

DO 68 Ks1,MZ 

AZ= 2。0*(DELTAT) *FM 

X(K)=(AZ)*X(K)/XENGY 

FREQ(K+1)=FREQ(K)+(1+0)/(AZ) 

TAU(K+1)=TAU(K)+DELTAT — 

CONTINUE 
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70 


80 


500 
101 
102 
103 


104 
105 
106 
107 
108 
109 
110 
111 
112 
113 
114 
115 


TABLE E-1 (Con't.) 


POWER SPECTRUM WRITE OUT INSTRUCTIONS 

WRITE (69103) 

WRITE (69107) XFACTsAZERO 

WRITE (69108) 

WRITE (69109) (TAUCK) 9А (К) s9FREQ(K) 9 X(K) sK=] ΜΖ) 

LABEL=4H 

ITITLE( 1)=8HJESBERG 

ITITLE(2)=BHPOWER SP 

ITITLE(3)=BHECTRAL D 

ITITLE(4)=8HENSITIES 

ITITLE(5)28H -- RUN 

READ (50115) ITITLE(6) 

00 80 J=7912 

ITITLE(J)=8H 
CALL DRAW(125»FREQ»X»0»0»LABEL»ITITLE»200».2»09.0»2929892»9»1»LAST) 
CONTINUE 

FORMAT (212,F84,55,2165I2) 

FORMAT(12F6.4) 

FORMAT (99H1 TUKEY SPECTRUM ESTIMATE S » 


10POWER SPECTRAL DENSTIY ANO AUTO-CORRe CALCS, 77) 


FORMAT (20H DATA CONTROL CARD //) 

FORMAT (30HO DATA RECORD-~-(P(1) #15100) //) 
FORMAT (10F 8-3) 

FORMAT (8H XFACT= F8e598H AlO)= Fl2.5 //) 
FORMAT (42H TAU(SEC) AN(TAU) FREQ(CPS) XN(FREQ) //) 
FORMAT (F9。36F10。5 ぁ 5F12。3。F10s5) 

FORMAT( 3X614) 

FORMAT (3X>016) 

FORMAT(31H1 ERROR IN SR DATA AT HEADING ) 
ЕОКМАТ( 40Х»13»7Х»016) 

FORMAT (22H  IDENT(N)»N*x1»NUMHDR  //) 
FORMAT ( 3XA8) 

END 
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824 
825 
828 
830 


822 
823 


826 


827 
829 


TABLE E-2 


SUBROUTINE DATALIDENT»MAX»KLIST»KFLAG) 
DIMENSION IBLOCK(1001)»KDATA(40009»2) 
COMMON KDATA» IBLOCK 
CALLING INSTRUCTIONS FOR SUBROUTINE DATA 
ARG le IDENT IS A 16 OCTAL DIGIT IDENTIFYING NUMBER 
SUPPLIED BY THE 160 PROGRAM WHEN TAPE WAS MADE. 
THE CALLING PROGRAM MUST SUPPLY THIS NUMBER IN ORDER FOR 
THE CORRECT BLOCK OF DATA TO BE RECALLED. | 
ARG 2„ МАХ IS THE NUMBER OF SAMPLES PER BLOCK AND MUST BE LESS: 
THAN 4000 DECe 
ARG 32 KLIST IS EITHER 1 OR 2 AND DECIDES WHICH SIDE OF 
KDATA YOU WANT THE DATA TO BE UNPACKED INTO 
ARG 4。 KFLAG IS AN ERROR FLAG RAISED BY DATAe IF =02NO ERROR 
LOGICAL UNIT 1 IS USED BY DATA TO FIND BLOCK. 
EACH TIME DATA IS CALLED IT FINDS THE DESIRED BLOCK ON THE TAPE 
AND UNPACKS IT INTO EITHER KDATA(Ms1) OR KDATA(M?2 ) 
IN THE CALLING PROGRAM DIMENSION KDATA(4000»2)» AND 
DECLARE KDATA COMMON, 
KDATA IS THE OUTPUT LIST AND IS REFERENCED BY KLIST 
SUBROUTINE DATA CALLS ON SR UNPACK AND SR FINDIT 
AFTER UNPACK THERE EXISTS ONE 1604 WORD/160 WORD 
ASSIGNMENT NOS USED IN DATA5,8225,823,82^,825,828,829 
FORMAT(48H1 ERROR OCCURRED IN FINDIT SUBROUTINE AT HEADER ) 
FORMAT(50X>016) ` 
FORMAT(33HO I HAVE UNPACKED DATA HEADED BY ,,016 ) 
ҒОКМАТ( 1Х%9184» 3Х» 016) 
CALL FINDIT(IDENTIMAX» IFLAG) 
IF(IFLAG) 822»823,822 
PRINT 824 $ PRINT 825»IDENT $KFLAG=1$ GO TO 829 
JMAX=MAX/4+1 
CALL UNPACK (IBLOCK» JMAX»*KLIST»KDATA(1IsKLIST)) 
KFLAG = O $ PRINT 828 » IDENT 
DO 826 M=1,4 | 
PRINT 830sM»KDATA(M»KLIST) 
J = MAX - 3 
00 827 М = Ј, МАХ 
PRINT 830eMsKDATA(MsKLIST) 
CONTINUE 
END 
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814 
815 
816 
818 
820 


805 
806 
807 
808 
811 
810 
812 
809 
813 
817 
819 


821 


TABLE E-3 


SUBROUTINE FINOIT( IDENT9»MAXs»IFLAG) 

DIMENSION KDATA(400092) sIBLOCK(1001) 

COMMON KDATAsIBLOCK 

FORMAT(53H1 A PARITY ERROR WAS DETECTED BUT RUN WAS NOT STOPPED) 

FORMAT(45H PARITY ERROR OCCURRED AT HEADER AS FOLLOWS ) 

ҒОКМАТ( 33Х% 016) | 

FORMAT(32H1 I HAVE LOCATED DATA HEADED BY »016 ) 

FORMAT (38H UNABLE TO LOCATE HEADING AS FOLLOWS ) 

MAX1=MAX/4 + 1 

ASSIGN 811 TO JUMP 

BUFFER IN (191)(IBLOCK(1)»IBLOCK(MAX1)) 

ІҒ(ММІТ»1) 806»807»808»810 | 

IF( IDENT—IBLOCK(1))805 98139805 

GO TO JUMP,(8115812) 

REWINO 1$ ASSIGN 812 TO JUMP $ GO TO 805 

IPAR=1 % GO TO 807 

REWIND 1$ GO TO 819 

PRINT 814% PRINT 8158 PRINT 816e9IDENTS GO TO 817 

IF(IPAR) 80998179809 

IFLAG = 0 $ PRINT 818»IDENT $ GO TO 821 

PRINT 820 

PRINT 8i6, IDENT 

IFLAG2] 

η NOS USED HERE»805»,807»,808,809,810,811,812,813,814,815 
8169817»818»819,820»821 

СОМТІММЕ 

END 
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UNPACK 


Ul 


AADRS 


41 


J2 
J3 

J4 
BADRS 


EXIT 


IDENT 
ENTRY 
SLJ 
SIU 
LIU 
LDA 
SAL 
ARS 
INA 
SAU 
LDA 
SAU 
LDA 
SAU 
INA 
SAL 
INA 
SAU 
INA 
SAL 
INI 
SIU 
SIL 
ENI 
ENI 
LDA 
LRS 
QRS 
STQ 
LRS 
QRS 
STQ 
LRS 
QRS 
STQ 
STA 
INI 
ISK 
SL J 
ENI 
ENI 
SLJ 


me Ree 
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N EP 


UNPACK 
UN PACK 
++ 
EXIT 
UNPACK 
0 

U1 

24 

-1 
AADRS 
xx 
BADRS 
1 

J4 

1 

J3 

1 

J2 

1 

J1 

2 
EXIT+1 
EXIT 
2 

0 

KH 

12 

36 

жж 

12 

36 

хх 

12 

36 

HH 

+ 

& 

жж 
AADRS 
xx 

жж 

κκ 


TABLE El} /, 
SUBROUTINE UN: CK 
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M 


ARGS IBLOCK»JMAX»KLIST+KDATA 
SUBR UNPACK CALLED BY SR DATA 
+ 


+ JMAX=MAX/4+1 。 MOST=JMAXー1 
GETS ADDR OF CALLING ARGS 

+U1 IS JMAX ADDR 

IBLOCK PACKED IN DATA 1001 wDS 
+ 

ADRS OF IBLOCK IN AADRS 
+FILLED IN BY 3 INSRUCTIONS BACK 
STORE JMAX IN LOOP COUNT 

+GO GET ADRS OF NEXT ARG 
+STORE ADRS OF KDATA(1sKLIST) 
ADVANCE ADRS 4 
+STORE ADRS OF KDATA(2sKLIST) 


+DITTO FOR KDATA(3sKLIST) 


+DITTO FOR KDATA(4,KLIST) 
SET CORRECT EXIT ADRS 
+STORE IN EXIT INST 

SAVE INDEX 2 

+FIRST DATA WORD IN IBLOCK 
FIRST WORD IN KDATA IS DATA 
+IBLOCK( J) 

SHIFT 4TH WORD INTO A REG 
*RIGHT JUST,SIGN EXTEND 
STORE IN KDATA (I+4»KLIST) 
+SHIFT 3RD WORD INTO Q 
RIGHT JUST »SIGN EXTEND 
+STORE IN KDATA(I+3»KLIST) 
SHIFT 2ND WORD INTO A REG 
+ RIGHT JUST» SIGN EXTEND 
STORE IN KDATA(I+29KLIST) 
+1ST WORD NOW RTeJUST +SIGN EXT 
1=1+4 

*ISK ON JMAX 

ο πο] REPEAT LOOP 

+RESTORE INDEX 1 

RESTORE INDEX 2 

+ JUMP OUT 


UNPACK 
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AADRS 


J1 


J2 
J3 
J4 
BADRS 


EXIT 


-EXECUTE. 


IDENT 
ENTRY 
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510 
110 
LDA 
SAL 
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INA 
SAU 
LDA 
SAU 
LDA 
SAU 
INA 


INA 
SAU 
INA 
SAL 
INI 
STU 
SID 
ENI 
ENI 
LDA 
LRS 
QRS 
STQ 
LRS 
QRS 
STQ 
LRS 
QRS 
STQ 
STA 
INI 
ISK 
SLJ 
ENI 
ENI 
SLJ 
END 
FINIS 
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UNPACK 
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UNPACK 
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TABLE E-4 (Con't.) 


ARGS IBLOCK,JMAX»KLIST»KDATA 
SUBR UNPACK CALLED BY SR DATA 
+ 


+ JMAX=MAX/4+1 y MOST=JMAX-1 
GETS ADDR OF CALLING ARGS 

+U1 IS JMAX ADDR 

IBLOCK PACKED IN DATA 1001 WDS 
+ 

ADRS OF IBLOCK IN AADRS 


+FILLED IN BY 3 INSRUCTIONS BACK 


STORE JMAX IN LOOP COUNT 

+60 GET ADRS OF NEXT ARG 
+STORE ADRS OF KDATA( 1 5KLIST) 
ADVANCE ADRS 

+STORE ADRS OF KDATA(2»KL IST) 


+DITTO FOR KDATA(3»KLIST) 


+DITTO FOR KDATA(4»KLIST) 
SET CORRECT EXIT ADRS 
+STORE IN EXIT INST 

SAVE INDEX 2 

+FIRST DATA WORD IN IBLOCK 
FIRST WORD IN KDATA IS DATA 
+IBLOCK(J) 

SHIFT 4TH WORD INTO A REG 
+RIGHT JUSTsSIGN EXTEND 
STORE IN KDATA (I+49+KLIST) 
+SHIFT 3RD WORD INTO Q 
RIGHT JUST »SIGN EXTEND 
+STORE IN KDATA(I+3»KLIST) 
SHIFT 2ND WORD INTO A REG 
+ RIGHT JUST» SIGN EXTEND 
STORE IN KDATA(I+2»KLIST) 
+1ST WORD NOW RTeJUST +SIGN EXT 
I=I+& 

+ISK ON JMAX 

J=J+1 REPEAT LOOP 

+RESTORE INDEX 1 

RESTORE INDEX 2 

+ JUMP OUT 
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TABLE E-5 
SAMPLE COMPUTER PROGRAM OUTPUT 
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TABLE E-5 (Con't.) 


SEE ERC T RUR Е 51 ШТА Т Е 5 


TULEV 


A(0)= © (9463 


. 99626 


XFACT= 


FREQ(CPS) XN(FREG) 


AN( TAU) 


TAU(SEC) 
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APPENDIX F 


DEFINITION OF AUTO-CORRELATION FUNCTION 
AND POWER SPECTRAL DENSITY 


It is well known that under certain conditions a periodic 


function of a real variable can be represented by a Fourier series: 


(=) = ມ (as COS arne + ba SIN раа 
Nz- o 2 T 
where T is the period of the signal and the coefficients are given 


by: 
= 


z 
E E n 
cu = fl) cos <= == + dé ; ba | 66) SIN at ce 
η 


-T 
t 


It follows from this that the power in a single periodic signal 
exists only at the frequency n/T cycles per second where n is any 
integer. The amplitude of the sinusoidal component of frequency 


n/T is a + n Thus the power in this component is 


n 
(a < + b. ^)/2. To express the power spectrum directly as a function 
of frequency, it is convenient to employ the Dirac delta function 


J (w-w,) which has the following formal properties: 


ὅ(ω-ω)- ο Fr Wk k, 


$ 5 (was) Ғ(ш) i f (w,) 


then the power spectrum can be expressed as: 
peur m 
p(w) = $ eb) (wem) 0 
N= -© 


The auto-correlation function of a periodic function is defined as: 


ly 


を (の) + + J eto) £ (€+ Ç ) dt 


າ 
2: 


Τι 


Computation of this quantity involves multiplication of the function 
by its time-shifted replica followed by integration over one period. 
Using the Fourier representation of equation (1) the correlation 


function of the periodic function f(t) is: 


ο o2 eo 


ວ | 
K Cc) "21 2 E cos am + Dn SIN 1 
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= = ST, 
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The Fourier transform of a cosine is given by: 


“Т. (со upt) = J cos шь ее Slow) (3) 


thus, it is readily seen from equations 1, 2 and 3 that for a 


periodic function 
- At 
φίω)- | R(te AT 


00): 5 | ре) еа 


- оэ 
that is, the power spectrum and auto-correlation function of a 
periodic signal constitute a transform pair. The equations are 


known as the Wiener-Khintchine relations. By letting =0 we obtain 


directly the Parseval theorem 


το e 
L^ eae c BOO d [ Ple)dw 
Л Y - oo 

This relation states that the power in a periodic signal can be 


obtained either by averaging the square of the signal over one period 


or by integrating the power spectrum [3]. 
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The definition of the auto-correlation function and power 
spectral density are quite clear for a periodic signal. However, 
a stationary random function is considered in this report. 

Let u(t) represent a random function of time. Normally a 
complete statistical description is necessary to define a random 
signal. For dynamic-load problems in aeroelasticity a time average 
concept may be used. The mean and the mean square of a function of 


time u(t) over a time interval 2T are: 


to tr 
$T 
— / coat 
а) = Z ) “* (€) dt 
ເ... 


If the time average is independent of to and T when T is 
sufficiently large, then u(t) is said to be stationary. 

Assume that the function u(t) is observed for a long period 
2T, and that it vanishes outside the interval 2T. Further assume 
that 2T is greater than any characteristic time of the random 
function. Hence up (t) represents all average statistical properties 


of the random function. 


time t 


ът 


Truncated Function 


19 


Since the truncation process assures absolute convergence and no 
finite discontinuities exist, the Fourier transform and inverse 


exist. Therefore, 


πω) = f « (9 e プ qe 


Urlt) = = fc (ເບ) ຂ イイ の 


since u(t) is real 


> 
Fy (w= “Fo ее) (5) 
* 
where 2 is the complex conjugate. 


The auto-correlation function is defined as: 
T 
c DIEI VISUS の M 
で (= 9 (G+ で ) = ジン ンー Zr | ult) u(er 7) at (6) 
-T 
substituting the transform definition of u(t+T ) 
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rearranging the order of terms, 
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since by setting T =0 


CAS 


eee ua 
46): 4 -# | oode = [pdo (5) 
- oo о 
By definition の id < Cos wT -a sm wt the imaginary 


part vanishes due to symmetry. 


4(t)- 2 f P (ὦ) cos wT edw 


By comparison with previous definitions 
co 
4 : + | (ГС) cos Mo 
— o 

It is quite evident that the power spectral density is a 
Fourier cosine transform of the auto-correlation function and for 
a random function this is quite significant since the auto- 
correlation function can be easily ວວ Аа and more clearly 
understood. 

Plotting ¢(w) against ™ produces a curve, the area under 
which represents the mean square of the random function as 
indicated in equation 8. The spectral character can appear either 
as a distributed form (type A) or as a dominant frequency (type B) 


or a combination of either. The figure below is in an approximate 


normalized form and illustrates the two types of responses. 
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